Incommensurate antiferromagnetic fluctuations in the two-dimensional Hubbard model 
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Commensurate and incommensurate antiferromagnetic fluctuations in the two-dimensional repulsive t - t'- 
Hubbard model are investigated using functional renormalization group equations. For a sufficient deviation 
from half filling we establish the existence of local incommensurate order below a pseudocritical temperature 
Tpc. Fluctuations not accounted for in the mean field approximation are important — they lower Tpc by a factor 
^2.5. 
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The two-dimensional Hubbard model (llSS] has attracted 
much interest in the past two decades because it is a candidate 
model for the Cu02 -planes in the high T^-cuprates and may 
exhibit d-wave superconducting order |4il5|] at finite chemical 
potential. The model shows other interesting order structures 
such as incommensurate antiferromagnetism which appears 
close to half filling. 

We focus on repulsive interactions U and not too large 
next- to-nearest neighbour hopping f , where the model is 
an antiferromagnet at half filling. Not so far away from 
half-filling a more complicated form of antiferromagnetism, 
namely incommensurate antiferromagnetism is suggested by 
mean field coniputations and numerical studies for finite sys- 
tems iSSillliSiniEHHIB]. Incommensurate 
antiferromagnetism is related to the existence of spiral mag- 
netic states which occur at large values of U 1 1611171118]. Ex- 
perimentally, incommensurate antiferromagnetism manifests 
itself in the peak structure of the magnetic structure factor 
which is accessible via neutron- scattering. It has been ob- 
served for a variety of high Tc^cuprates, for experimental and 
numerical results see CIEMIHBIMIiI. 

In the temperature region where local incommensurate an- 
tiferromagnetic order supposedly sets in, the eflTective inter- 
action between the electrons is large such that perturbative 
methods are not reliable. Collective fluctuations of electron- 
hole pairs in the antiferromagnetic channel play an important 
role. Since they are omitted in a mean field treatment one 
may doubt whether the mean field results for incommensurate 
antiferromagnetism are reliable. For these reasons we inves- 
tigate the issue of incommensurate antiferromagnetism by a 
method that is intrinsically non-perturbative and includes ef- 
fective collective bosonic fluctuations, namely the functional 
renormalization group for the "flowing action" (or "average 
action") fldjll]. For this scale dependent eflfective action 
(or coarse grained free energy) the scale k indicates an infrared 
cutoflT such that only fluctuations with momenta larger than k 
are eflTectively included. (Finally, one is interested in the limit 
^ ^ 0, where Fy^^o equals the eflfective action — the generating 
functional of IPI-correlation functions — including all fluctu- 
ations.) We will work in a version where the dominant col- 
lective bosonic fluctuations are represented by bosonic fields 
II28II29'] . Our model is equivalent to the purely fermionic Hub- 
bard model from which it is derived by means of a Hubbard- 
Stratonovich transformation (SOl |3l|] . Earlier studies employ- 
ing the present framework have focused on the temperature 



dependence of commensurate antiferromagnetic order ||28 
the Kosterlitz-Thouless transition in a more general class of 
Hubbard-type models |32], and the generation of a coupling 
in the J- wave superconducting channel |33]. The role of in- 
commensurate antiferromagnetic fluctuations was not taken 
into account in this earlier work. Functional renormalization 
group treatments of the Hubbard model are more often given 
in a purely fermionic formulation, see L34l[35l[36[[37l 38,.39i]. 
Ref. m is of particular interest since, in accordance with 
the results described here, it also reports on a region in the 
phase diagram where incommensurate spin density fluctua- 
tions dominate. The present paper also includes an indepen- 
dent computation of the size of the incommensurability that 
occurs. 

Our ansatz for the flowing action includes contributions for 
the electrons, for the bosons in both the antiferromagnetic 
and d-wave superconducting channels, and for interactions be- 
tween fermions and bosons: 



The collective field ;^ = (a, J, J*, </^*) describes fermion 
fields the "antiferromagnetic boson field" a and the 

complex field d a finite expectation value of which signals 
J- wave superconductivity. The fermionic kinetic term 



involves the inverse fermion propagator 



(2) 



(3) 



where ^(q) = -i^-2t(cosq:f+cosqy)-4f cos qx cos qy depends 
on the chemical potential fi and the nearest and next-to-nearest 
neighbor hopping parameters t and f of the Hubbard model. 
We employ a compact notation X = (t, x), Q = {m, q), 



5(X-X') = 5(r-r'Mx,x', 
S(Q - Q') = J36„,„,i2nf6'-^\q - q') . 



<1_ 



(4) 



where all components of X or 2 are measured in units of 
the lattice distance a or a"^. The discreteness of the lattice 
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is reflected by the 27r-periodicity of the momenta q. A scale 
dependent fermionic wave function renormaUzation is in- 
cluded in Eq. ©. 

The purely antiferromagnetic bosonic term is described by 
a kinetic term and a local eff'ective potential 

Q X 
where we employ a quartic efl'ective potential Ua for a: 

Ud^] = mla+ ^Aaa^, (6) 

with a = a^/2. The kinetic term Pa involves the Q-dependent 
part of the inverse antiferromagnetic propagator and there- 
fore contains the essential information about difl'erent kinds of 
magnetism. Our treatment of this term is discussed in detail 
below. Local antiferromagnetic order in domains of size 
is signalled by a minimum of F^,;^ for a(2) ^ 0. For 2 = this 
describes commensurate antiferromagnetism, while a mini- 
mum for non- vanishing q in g = (0, q) indicates incommen- 
surate antiferromagnetism. A Yukawa-like interaction term 
couples the bosonic field to the fermions, 

rF^,k = -hai:K,Q,Q' S(K^U-Q^ Q) (7) 

where the momentum vector n is given by H = (0, tt, tt). 

The bosonic field d is associated to Cooper-pairs in the 
d-wave channel. It is described in more detail in (H, [sH] 
where the exact form of Yp^ can be found. In this note 
we include the eff'ect of d-wave fluctuations on the flow of 
the fermionic and "antiferromagnetic" part of F^. We use 
'^d,k =]:Qd\Q)PAQ)d{Q) + U^Md""] with J*] = 
m^S-\-AdS^/2 where S = d*d. Here we focus exclusively on the 
emergence of (either commensurate or incommensurate) mag- 
netic order which occurs in the vicinity of half- filling. The 
emergence of J- wave superconducting order at larger values 
of \fi\ will be discussed in detail in a future publication. No 
superconductivity has been detected in the region of the phase 
diagram studied here. 

The dependence of the flowing action on the scale k is de- 
scribed by an exact flow equation [26]. Our ansatz ([T]) approx- 
imates the solutions to this functional difl'erential equation. At 
the microscopic scale ^ = A the flowing action must be equiv- 
alent to the microscopic action of the Hubbard model. Since 
we want to eliminate the (constant) four-fermion coupling 
U at k = A which, of course, has no contributions exhibit- 
ing J-wave symmetry, the repulsive interaction between the 
fermions must be contained in the antiferromagnetic Yukawa 
coupling ha. In the bosonized picture one has, instead of the 
original four-fermion coupling U a boson-mediated interac- 
tion term hl/fhl which must be chosen proportional to U. 
Since an additional sum over spin directions has to be car- 
ried out it has to be chosen as U/3. Thereby we have simply 
transcribed the original model into an equivalent one using 
bosonic language. Since the original model does not contain 
any eight electron terms, no quartic bosonic coupling Aa can 




FIG. 1: Mean field kinetic term Pa(0, q)/t of the a-boson as a func- 
tion of space-like momenta for U/t = 3 and f = 0. In Figure (a) 
= and T/t = 0.205, in Figure (b) ji/t = -0.27 and T/t = 0.0435. 
Both temperatures are mean field critical temperatures. 

exist at the UV scale ^ = A. In sum, a set of possible "initial 
conditions" for the flow of the coupling constants is given by 

fhlW = u/3 , haW = U/3 , ^Ia = , Pa(Q)\A = 0, 

rFd\A = 0, F^|a = JV, Z^|a = 1. (8) 

These values specify the action Fa (or, equivalently, the 
Hamiltonian) at the microscopic level. 

At the microscopic scale A, the action for a is Gaussian 
such that a can be "integrated out" by solving its field equation 
as a functional of i//. The J-boson decouples and becomes 
irrelevant. This demonstrates that Fa indeed coincides with 
the purely fermionic Hubbard model with repulsive coupling 
U>0. 

We still have to specify the truncation for the kinetic term of 
the a-boson Pa(Q). This is a central object of this paper, since 
incommensurate antiferromagnetic fluctuations will dominate 
if the minimum of Pa(0, q) occurs for nonzero q. In order to 
gain some first information about the general shape of Pa we 
compute the mean field contribution from the fermionic loop 

hi 

^Pa(Q) = > + (G ^ -2) . (9) 

The general features of this mean field contribution are used 
in order to motivate the form of the bosonic propagator in our 
truncation. We observe close to half filling two qualitatively 
diff'erent situations. At half filling and for sufficiently high 
temperatures also close to half filling there is a pronounced 
minimum at q = 0, see Fig. [I](a). However, away from half 
filling the picture is diff'erent for sufficiently low temperatures, 
see Fig. [l](b). In the center at q = there is a local maximum 
and there are four minima at positions 

qi,2 = (±^,0), q3,4 = (0,±^), (10) 

where ^ is a function of T, yu, and f. This is a manifestation 
of the dominance of incommensurate antiferromagnetic fluc- 
tuations. 

Once the minimal value of the inverse bosonic propaga- 
tor Lp^(0, q)-Fm^ at zero frequency becomes smaller than 
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zero, the minimum of the free energy can no longer occur 
for (a(2)) = 0. One rather expects spontaneous symmetry 
breaking with a non-zero expectation value of (a). As long 
as the minimum of ^^(0, q) is located at q = 0, the order 
parameter (|a|) ~ ^(q) indicates commensurate antiferromag- 
netism. However, for a minimum at q = q^ 9^ the in- 
commensurate antiferromagnetic order breaks further lattice 
symmetries. One of the pairs of minima (fTOl) is selected and 
the symmetry of rotations by njl around q = in momen- 
tum space is spontaneously broken. The spins change sign 
between neighboring lattice sites only in one direction, the 
x-direction say, whereas in the orthogonal direction the peri- 
odicity corresponds to some momentum 7i±q. The state with 
(a) = becomes unstable when (Pa)min = -^l- In case of 
a second order phase transition this occurs for the mean field 
critical temperature T = Tmfc- Figures [T] (a) and (b) corre- 
spond to mean field critical temperatures. Note that the sys- 
tem selects one of the pairs qi 2 or q3,4 since a(X) is a real 
field. Therefore the system remains symmetric with respect to 
reflection about the axes. 

We are interested in whether incommensurate antiferro- 
magnetism persists if bosonic fluctuations are included. Tak- 
ing into account bosonic fluctuations, the critical tempera- 
ture vanishes in the infinite volume limit due to the Mermin- 
Wagner theorem. The destruction of local order by the long 
range fluctuations of the Goldstone bosons (antiferromagnetic 
spin waves) is only a logarithmic efl'ect, however. For a probe 
of finite macroscopic size antiferromagnetic order can be ob- 
served and the critical temperature is nonzero |29]. Here the 
eff'ective critical temperature Tc is defined such that for T < Tc 
the typical size of ordered domains exceeds the macroscopic 
size of the probe /. In other words, (a(^)) diff'ers from zero for 
kph = if T < Tc, while for T > Tc one has (a(^)) = 0. 

In this note we only study the pseudocritical temperature 
Tpc which marks the onset of local ordering corresponding to 
a minimum of the flowing action Tk for a(0, q). Above this 
temperature, (|a|) = holds on all scales of the renormaliza- 
tion flow. For T < Tpc local order sets in for ^ = kc > 0. 
In case of incommensurate antiferromagnetism we expect the 
formation of domain walls between regions where q points in 
the X- or j-direction. This constrasts with commensurate anti- 
ferromagnetism where only a continuous symmetry is broken 
for a 9^ 0. For k < kc the flow should then be continued in 
a regime with nonzero a in order to account properly for the 
Goldstone boson fluctuations. This has been investigated for 
the commensurate case in 12?], but is not yet implemented for 
the incommensurate case in the present note. We note that Tpc 
is the equivalent of the mean field critical temperature. For 
Tc < T < Tpc the electron propagator does not exhibit a true 
gap, but it is suppressed for momenta corresponding to the 
inverse of length scales for which local order is present. 

Inspired by the shape of Pa in the mean field approximation 
we approximate the kinetic term for the antiferromagnetic bo- 
son by 



(11) 




FIG. 2: Renormalization flow of the gradient coefficient A^^t for 
U/t = 3 and fi/t = -0.12 according to Eq. ([Tsl l in the parameter 
regime where we have commensurate antiferromagnetism and there- 
fore dkq = 0. The solid line corresponds io T jt - 0.08, the long- 
dashed line to Tjt = 0.07, and the short-dashed line to Tjt = 0.058. 



choose for F{q) 



Z)^[q]^ 
Z)2 -h [qp 



(12) 



Here [q]^ is defined as [q]^ = q\ + for qi e [-71,71] and 
continued periodically otherwise. For small q^ the quadratic 
approximation Pa = A^q^ describes a linear dispersion rela- 
tion for the composite bosonic field, oj = ^fAJZ^\q\, while 
for q near the boundary of the Brillouin zone the momentum 
dependence of Pa 'levels off'' as in Figs. [I](a), (b). For a suit- 
able choice of Aa and D the shape of the mean field result for 
Pa is well reproduced. Of course, due to the important contri- 
butions of bosonic ffuctuations beyond mean field theory, the 
actual values of Aa and D will diff'er substantially from the 
mean field values. 

Within the functional renormalization group approach, we 
describe the scale dependence of the bosonic kinetic term by 
ffow equations for the parameters Aa and D. For these pur- 
poses we define the gradient coefficient Aa by 



1 , 
A. = -^P.(0,/,0)^ 



(13) 



with § = in the commensurate case. The shape coefficient 
D is computed as 



1 



(Pa(0,7r,7i)-Pa(0,q,0)). 



(14) 



where for the case of commensurate antiferromagnetism we 



The ffow equations for Aa and D can be extracted by inserting 
our truncation in the exact ffow equations for the kinetic term 
CD). 

During the renormalization ffow the gradient coefficient Aa 
ffrst increases, starting from = at the scale A. At half 
ff lling and in the proximity of half ff lling for sufficiently high 
temperatures, Aa either increases monotonically or at least re- 
mains larger than zero on all scales k < A, see Fig. O The 
minimum of Pa occurs for q = and commensurate antifer- 
romagnetic ffuctuations dominate. 



FIG. 3: In (a) the mean field approximation for bosonic kinetic term 
Pa(0,qjc,0)/t is shown as a function of the x-component of spatial 
momenta. Parameters are U/t - 3, ju/t - -0.35, f = and T/t = 
0. 1. Fig. (b) shows the same quantity according to our approximation 
given by Eqs. iTTl and iT6i with the values of A^, q and D drawn 
from the mean field computation. 



However, for low enough temperatures and at sufficient dis- 
tance from half-filling, Aa becomes zero on a certain scale. If 
we continued to evaluate Aa for q = it would decrease to 
negative values for lower scales. This situation corresponds to 
the case of incommensurate antiferromagnetism. The ansatz 
for the function F(q) given in Eq. ([T2l) is no longer suitable. 
One has to allow for the existence of minima of Pa,k(^, q) at 
nonzero 0. The ansatz ([TT1) for the inverse bosonic prop- 
agator employs now for F(q) 



D^^F(q,q) 



(15) 



The quadratic momentum dependence of the numerator in 
([T2I) is replaced by an expression which is quartic in momen- 
tum and explicitly includes the incommensurability q: 



[q]¥ + 4[^j'[^/). 



(16) 



The first term in F vanishes for [q]^ = q^ and suppresses the 
propagator for [q]^ 9^ q^. The second term favours the minima 
([TOl) as compared to a situation where rotation- symmetry in 
the qx - ^3; -plane is preserved. The prefactor is determined by 
Eq. ([T3]) . For q ^ one has Aa ~ q^ such that Pa becomes 
quartic in q. We compare in Fig. [3] the kinetic term Pa(0, q) 
in mean field theory with the approximation from our ansatz 
which shows satisfactory agreement. 

In Fig. m a typical flow for Aa and q in the incommensu- 
rate regime is displayed. For scales below the scale where Aa 
becomes zero, q increases to a finite value and Pa(0, q) has 
four degenerate minima at positions given by Eq. ([TOl) . The 
solution q of Eq. (|2Q1) at the end of the flow corresponds to the 
position of the minimum, e.g., at the positive ^^-axis .We next 
specify the flow in more detail. 

The regulator function R^iQ) for the antiferromagnetic fluc- 
tuations should be adapted in order to allow for the domi- 
nance of incommensurate antiferromagnetism. We employ, 
similarly for the commensurate and incommensurate case, 

Rt(Q) = • (k^ - FcM^ q))e(k^ - FcM^ q)) , (17) 



respectively. This generalizes the cutoflT chosen in [.i33|] . 
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FIG. 4: Renormalization flow of the gradient coefficient Aa,k and the 
incommensurability q according to Eqs. fTSl ) and (l20l ) dXfilt - -0.12 
and T It - 0.05, where incommensurate antiferromagnetism domi- 
nates. The solid line shows Aa decreasing to zero dX - \nklt - 3.09 
and increasing again when the incommensurability q (short-dashed 
line) sets in. 



The flow equation for the gradient coefficient is obtained by 
taking appropriate derivatives in one of the minima 



dkAa 



Q 



1 



dPp\{Q)P^p{K+Q + IV) 
1 



dp P%{Q)P%{K+Q + IV) 



, (18) 



Pf(Q) + 
The first 



where K = (0, /, 0). We have defined PfM 
R^(Q), with fermion cutoff R^ chosen as in 
term in (fTSl) results from the change of the infrared cutoflf in 
the fluctuations. The symbol dk means a formal derivative 
with respect to the cutoflf function R^. The second term in 
(fTSl) reflects the shift of the location of the minimum of Pa 
at (q, 0) and is absent if commensurate fluctuations dominate, 
q = 0. 

A flow equation for the position of the minima q is derived 
from the condition 



dqx 



(19) 



Taking the scale derivative of this equation one obtains the 
flow equation: 

(«^/-.,.(o.„u„,4L^.,.(o,„u,, ,20, 



dk\qdqx 



= 0. 



Flow equations for the other running couplings 
Zf,m^,Aa,ha,m^,Ad,hd,D are not given explicitly here, 
see 133]. 

We now turn to the results obtained in our renormalization 
group scheme. An overview of the occurence of incommen- 
surate antiferomagnetism is given in Fig. [S] showing pseud- 
ocritical temperatures Tpc for the diflTerent kinds of antiferro- 
magnetic order in the presence of vanishing (upper panel) and 
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FIG. 6: Renormalization group results for the incommensurability q 
as a function of T for U It - 3, f - and jilt - -0.12 (solid line) 
and jdlt = -0.105 (dashed line). 



0.08 



0.06- 

Tlt ' 

0.04 




0.00' ' ■ ' ■ ' ' 

-0.3 -0.2 -0.1 0.0 0.1 

/jIt 

FIG. 5: Renormalization group results for the pseudocritical tem- 
perature Tpcit as a function of It, given by the solid (commensu- 
rate) and dashed (incommensurate) lines. Results are displayed for 
Ult -?),t' =0 (upper panel) and t' It - -0.05 (lower panel). 



nonvanishing (lower panel) next-to-nearest neighbor hopping 
f . The solid line signals the onset of local commensurate, the 
long-dashed line the onset of local incommensurate antifer- 
romagnetic order. Below the short-dashed line there is no lo- 
cal magnetic order but incommensurate fluctuations dominate. 
Below the point where the short-dashed line terminates at low 
temperatures, numerical solutions to the flow equations, as we 
have implemented them numerically, are no longer reliable. 
For both vanishing and non- vanishing t\ one observes com- 
mensurate antiferromagnetism for a certain range of chemi- 
cal potential yu, while for smaller and larger values of ji in- 
commensurate fluctuations begin to dominate. For finite t\ 
however, the pseudocritical curve is no longer the same for 
positive and negative fi but, for negative t\ is shifted to more 
negative values of jd. 

The pseudocritical temperature is found to be substantially 
lower than according to the mean field computation. For U = 
3t, f = and yu = 0, for example, the mean-field computation 
gives TuFcIt = 0.205, while we find Tpc/t = 0.0745 when 
one takes into account bosonic fluctuations. By reducing the 
interaction, the shape of the pseudocritical curve remains the 
same but local order emerges only at lower temperatures. 

With decreasing temperature the tendency towards incom- 
mensurate fluctuations is increased, which can be demon- 
strated by studying the dependence of ^ on T at fixed chem- 
ical potential. It is shown for ju/t = -0.105 and j^/t = -0.12 



in Fig. [6l For large enough temperatures one has q = 0, while 
below some yu-dependent temperature incommensurate anti- 
ferromagnetism sets in. The temperature where this happens 
is indicated by the short-dashed line in Fig. [5] (upper panel). 
For smaller T the value of q increases, the final point of the 
jj/t = -0.105-curve at low temperature corresponds to the 
long-dashed line in Fig. \5\ 

As one can see from the curve representing fi/t = -0.12 in 
Fig. [6l at small temperatures the size of the incommensurabil- 
ity is approximately constant. Therefore we compare our re- 
sult to the zero-temperature result obtained by saying that 
q = 2 arcsin(|yu|/20 (which has also been used in the fermionic 
RG computation given in |36D). For ji/t = -0.12 this for- 
mula gives q ^ 0.120 whereas we find q ^ 0.132. By taking 
into account fluctuations the incommensurability seems to be 
slightly enhanced. Agreement with the results displayed in 
1 12] obtained by means of the composite operator method is 
also satisfactory. 

A dominance of incommensurate antiferromagnetic fluctu- 
ations can be observed in the momentum dependence of the 
magnetic susceptibility and the bosonic occupation number. 
The susceptibility is given by the bosonic propagator at zero 
frequency Pa^(0, q), while the occupation number is obtained 
by an additional sum over bosonic Matsubara frequencies, 
Haiq) = TYjajB(Pa((^B,q))~^' Fig. [7] shows that for param- 
eters where the bosonic mass is small, here fhl/U ^ 10"^, and 
thus close to the onset of local incommensurate order, both the 
magnetic susceptibility and the bosonic occupation number 
are peaked at qx = ±q, signalling that incommensurate fluc- 
tuations strongly dominate. The situation is completely anal- 
ogous for the -dependence of the susceptibility at qx = 0, 
whereas both quantities do not have such a pronounced peak 
structure along the Brillouin zone diagonal. 

In those regions of the phase diagram in which (either com- 
mensurate or incommensurate) antiferromagnetic order exists 
on a certain legth scale k our truncation becomes inapplica- 
ble in the regime below k. The simplest way of obtaining a 
glimpse at these regimes is by means of a mean field analysis, 
so before closing the discussion we briefly address this prob- 
lem. A more extensive mean field treatment, if only with re- 
gards to the commensurate case but including a nonzero next- 
to-nearest neighbor hopping t\ is given in 1,40^ . Here one has 
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FIG. 7: Fig. (a) shows the spin susceptibility (Pa(0,qx,0)/t) and 
Fig. (b) the bosonic occupation number na(qx,0) for fi/t = -0.12, 
T/t = 0.01, U/t = 3 and f = according to the renormalization 
group computations. Both curves are given as a function of spa- 
tial momentum in x-direction. A peak ai =0 would signal the 
dominance of commensurate antiferromagnetism. The actual peaks, 
located at § = +0.132, indicate incommensurate antiferromagnetism. 



to take into account that the periodicity of a system in the Neel 
state is changed resulting in a new "magnetic" Brillouin zone 
whose boundaries are given by the lines between the (±7r, 0) 
and (0, ±7i) points. Correspondingly, the mean field disper- 
sion relation for a nonzero gap parameter A = ha{\sL\} has two 
branches 



^±(P) = ^|^(p)+^(p + 7r): 



a/(^(p)-^(P + ^))' + 4A2|21) 



which, for finite f, lead to an interestingly structured eff'ective 
Fermi surface enclosing hole pockets around (±7r/2, ±7i/2) 
and electron pockets around (±;r, 0) and (0, ±7r), see the ex- 
ample drawn in Fig. [8] (a), for further details see e. g. ||4QIi . 






J 


f 







(b) 



FIG. 8: Mean field eff'ective Fermi surfaces for \ilt - -0.6, f It - 
-0.2 and gap parameter Ajt - 0.1. Fig. (a) shows the commensu- 
rate case q - ^ where the Fermi surface exhibits hole- and particle 
pockets at the magnetic Brillouin zone boundary. In Fig. (b), the 
remainders of the eff'ective Fermi surface are shown for a nonzero 
incommensurability q - 0.3 along the x-axis. 



In the presence of a nonzero expectation value <a(q)) with 
q 9^ 0, i. e. in the presence of incommensurate order, the in- 
verse of the fermionic mean field propagator at zero frequency 
has contributions from Eqs. ©(with Zy - \) and (|7]) and is 
given by 



X'O- 



= ^(q)^(q-qO (22) 
(^(q - q' - TT + q) + ^(q - q' - TT - q)) 
q3,4 as defined in Eq. (fTOb . The analogue 



V2^ 

with q = {{12 or q : 
of the Fermi surface corresponds to the zero eigenvalues of 
Pp. However, the corresponding eigenmodes are no longer 
momentum eigenstates. Nevertheless, if the gap parameter 
A = |A| is nonzero but small, many eigenvalues of Pfi^i, qO 
have most of their support each at a single momentum p. This 
concerns all those momenta p for which the condition 



A « |^(p -h TT + q)|, |^(p + TT - q)| 
is fulfilled. With respect to these momenta the equation 
1 1 



2 \^(p-h7r + q) ^(p + TT-q) 



= 



(23) 



(24) 



defines an eff'ective Fermi surface which is obtained by (ap- 
proximately) diagonalizing Pfi^, qO for small A. For large 
enough A the eff'ective Fermi surface vanishes completely be- 
cause the number of solutions to Eq. (l24l) that satisfy the con- 
dition (|23]) rapidly goes down. In Fig. [8](b) the eff'ective Fermi 
surface is shown for the incommensurate case with an order 
parameter <a(q)) where q = qi,2, i. e. the incommensurabil- 
ity is along the x-axis. The symmetry of rotations by tt/ 2 is 
manifestly broken. 

To summarize, we have shown that incommensurate anti- 
ferromagnetic order in the two-dimensional Hubbard model 
persists if bosonic fluctuations are taken into account. This 
phenomenon occurs at least in the form of local order for tem- 
peratures smaller than the pseudocritical temperature shown 
in Fig. \5\ We speculate that for T ^ the size of the incom- 
mensurate domains grows beyond the size of typical macro- 
scopic probes, but this remains to be shown. If magnetic fluc- 
tuations play a role in the generation of d-wave superconduct- 
ing order, the efl'ect of incommensurability has to be taken into 
account. 
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